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The Wigner function, which provides a phase-space description of quantum systems, has various 
applications in quantum mechanics, quantum kinetic theory, quantum optics, radiation transport 
and others. The concept of Wigner function has been extended to quantum fields, scalar and 
electromagnetic. Then, one deals with the Wigner functional which gives a distribution of field and 
its conjugate momentum. We introduce here the Wigner functional of fermionic fields of the values 
in a Grassmann algebra. Properties of the functional are discussed and its equation of motion, which 
is of the Liouville's form, is derived. 

I. INTRODUCTION 

A description of a quantum mechanical system is usually formulated in a position or momentum space. A prob- 
abilistic description in a phase space spanned by positions and momenta is not possible because of the uncertainty 
principle. The quasiprobabilistic approach can be achieved by using the Wigner function [![, which for a system of 
single degree of freedom, is usually defined nowadays as 



W(x,p;t) d ^ 



= dut 



-ipu 



x + u/2\p(t)\x-u/2), 



(1) 



where p(t) is the time-dependent density operator and \x) is the position eigenstate. The Wigner function is real but 
it is not everywhere positive and consequently the approach is not probabilistic but only quasiprobabilistic. Various 
quantum mechanical problems can be formulated in terms of the Wigner function Q. It is particularly useful to 
discuss the classical limit of quantum mechanics. In case of on many-body systems, the Wigner function plays a role 
of quantum analog of classical distribution function, see e.g. [3j], and it is thus the key object of quantum kinetic 
theory. The phase-space methods are particularly well developed in the theory of radiation transport and quantum 
optics, see the monographs [J, |5j. 

The concept of Wigner function was extended to quantum fields in Q where a real scalar field was discussed. The 
Wigner functional is defined in a full analogy to Eq. ([1} that is 



W[$,IL;t] 



del' 



X>y(x) exp 



(fx n(x)<^(x) 



($(x) + ^(x)|p(0|$(x)-^(x)} ! 



(2) 



where <I?(x) is the time independent scalar field in the Schrodinger picture, x £ I 3 is the coordinate variable and 
/ Vip(x) . . . denotes the functional integral. The Wigner functional was shown to be a density in an infinitely 
dimensional phase space spanned by ^(x) and its conjugate momentum II(x). Properties of II; t] were discussed 
and its equation of motion was derived. The scalar fields in thermal equilibrium and a system undergoing a phase 
transition due to the falling temperature were also studied. Later on the Wigner functional of electromagnetic field 
was analyzed 0] with the thermal, coherent and squeezed states being discussed. 

Our aim here is to introduce the Wigner functional of fermionic fields. Since the fields, which take their values 
in a Grassmann algebra, are very different from classical (bosonic) fields, it is a priori not obvious that the Wigner 
functional of fermionic fields, which has a meaning of phase-space density, can be constructed. We show, however, 
that, in spite of very different mathematical structure, the Wigner functional of fermionic fields can be defined as 
that of bosonic ones and that the two appear to be rather similar. Even more, at a superficial level the corresponding 
formulas look almost the same and the differences are merely in the mathematical background. It makes us to 
adopt a specific style of presentation with many - hopefully not too many - technical details included, as it is rather 
difficult to expose the differences and similarities of Wigner functional of bosonic and fermionic fields not referring to 
technicalities. 

In the subsequent sections the Wigner functional is defined and its properties as a phase-space density are discussed. 
The equation of motion, which has a from of the Liouville's equation, is also derived and its solutions are given. Some 
remarks close the paper. 
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II. DEFINITION OF WIGNER FUNCTIONAL 

Let us define the Wigner functional of fermionic field in the same way as in the case of bosonic field @ that is 

W[%U;t] = J 2ty(x) cxp[-* J d 3 x IT(xMx)] <*(x) + ^(x)|p(t)|*(x) - ^(x)), (3) 

where \&(x) is the time independent fermionic field in the Schrodinger picture and J Vip(x) . . . denotes here the 
functional Berezin's integral which is extensively discussed by its founder in Q . 

To check whether the definition makes any sense we consider for the beginning the Gaussian functional with 

(*l(x)|p|* 2 (x)} = e -H d3a; [*i( x )^*i( x )+*2(x)A* 2 (x)^ (4) 

where VPi, ^2 are the two-component Weyl spinors and A is an antisymmetric matrix (A T = —A). If the matrix A is 
not antisymmetric, only its antisymmetric part contributes to @, as squares of fermionic fields vanish. The formula 
(J4j) substituted into the definition ((3]) provides 

To compute the functional integral over <p(x), the field is discretized y(x) — > ipi, ip%, . . . tp n and using the formula 

dx m ■ .-dxi#m ■ ..dip 1 e XiM ^ +hiXi+Pi ^ = detAf e'^' 1 ^ , (6) 



which is not difficult to prove following the methods described in Q , the Gaussian Wigner functional is calculated as 

W[#,II] = C e -/ rf3 ^[*( x )^*( x )-n(x)A- 1 n(x)] j ^ 

where C is a constant. The result ([7]), which is rather expected as fully analogous to the case of bosonic fields, suggests 
that the definition @ makes sense and is worth further consideration. 

III. PROPERTIES OF WIGNER FUNCTIONAL 

The Wigner functional ([3]) is expected to represent a density in a phase space spanned by *f> and II. Then, there 
should hold the relation 

(o(*,n)) = (o(*,n)), (8) 

where 0(^,11) is a function of \E' and II. We use the hats to distinguish the field operators from their 'classical' 
counterparts. The word classical is written in the quotation marks as the anticommuting fields are not classical in 
the usual sense. The average of the field operators is defined as 

<0($, ft)) d = f Tr[p(t)0($, ft)] (9) 

and that of 'classical' fields in the following way 

(0(*,n)) d ^ f i Jv^^O^,IS)W[%U;t] (10) 

with 

f £>n r 

Z = J — W[^,U;t}. (11) 
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A. Computation of Z 



We start with the computation of Z presenting here some technical details which will be mostly skipped in the 
subsequent sections. Substituting the Wigner functional ([3]) into Eq. (jTTJ) and performing the discretization of the 
fields, one has 



, dll n _i dUi 

Z = I d^nd^n-x . . . d^i— . . . ——dip n dip n ^i . . . dipt exp 

Ztt Ztt Ztt 



fc=i 



x<*(x) + ^(x)|p(t)|tf(x)-^>(x)}, 



(12) 



where = ^(x.;), IT = ll(xi) and <fi = y(xi) with i = 1, 2 ... n; AV is the volume of elementary cell coming from 
discretization of R 3 . The equality holds in the continuum limit which is implicitly assumed here. Now we are going 
to take the integrals over IT which requires moving dUi right of d(fi. If n is even, the number of interchanges of 
Grassmann elements is also even and the sign of the whole expression is not altered. When n is an odd number, the 
sign changes. One computes the integrals over IT expanding the exponential and taking into account the n— th term 
which is the only one contributing to the integral. Then, one obtains 



dH n dTl n -i dlli 



2?r 2tt 



'Ztt 



exp 



Yk 



k=l 



iAV) n fdUndU^ dlli r^A 



-[En 

2tt 2tt 2tt 



kfk 



k=l 



= (-iAV) — — IIi^i H 2 (p 2 ■■■lL n (p ri 

Ztt Ztt Ztt 



(13) 



where the factor (—1)" , which occurs due to interchanging cflT and ipj, is included and the function 

5 (n) (<pi,^2, = <Pl<P2 ■ ■ -fn, (14) 

appears to be the Dirac delta-like function in the n dimensional Grassmann algebra. Indeed, 

d(f n dipn-l ■ ■ ■ dif! S^(ipi,ip2, ■ ■ ■ (fin) f(<Pl,<f2, ■ ■ ■ fn) = fa, (15) 

where fk with k — 0, 1, 2, ... n is defined through the decomposition of an arbitrary function f((pi, <fi2, ■ ■ ■ fn) as 



f(<Pl,<P2,---<fn) = fo + 22ft(i)<Pi+ 22 f%{h,i2) VhVia H h ^ fn{h,h, 



Jn) l Pii ( Pi2--- ( Pi n - (16) 



i=l 



The continuum limit of S^ n '((pi, f2, ■ ■ ■ fn) is denoted as 5[ip] and we have the identity 



'Ztt 



exp 



i / d xU(p 



C s 8[tp], 



(17) 



which is repeatedly used further on. It is a complex mathematical issue to determine the constant C$ which is, 
however, irrelevant for our considerations. It will be shown that the constants of this type drop out when the final 
average (ITUl) is computed. 

Using the result (|T3"|). Eq. ([12) provides 



Z = C J Z>tf(#|p(i)|*) = CTr[p(t)], 

Tr[p(t)}= J P*W(*)I*>- 
Now we are going to prove the relation (jSJ considering step by step three special cases. 



(18) 



(19) 
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B. C(*,f[) depends only on * 



We first consider 0(^,11} = *(x). Since squares and higher powers of fermionic fields vanish, there is no other 
local operator made of \I/(x) except a linear function of ^(x). Therefore, we compute 

1 f 1 



Tr[p] 



P*(* 1/5*1*) 



Tr[p] 



P* * (*|/3|*) 



(20) 



The last equality holds because the states |*) are by definition the eigenstates of * that is *|*) = *|*). Computing 
(*) exactly as Z, one finds 



P** (W|/5|#) = (*). 



(21) 



As seen the constant C cancels out and consequently we will ignore constants of this type further on. If one considers 
an average of a non-local product , I'(x 1 )*(x2) . . . f(x„), the result obviously is 



(*(xi)*(x 2 ) . . . $(x„)) = (*( Xl )*(x 2 ) . . . *(x n )). 



(22) 



C. e>(*,n) depends only on II 



As previously, we first consider C(*,il) = li(x) 

1 



^n(n|p|n), 



where 



Tr[p] 



PIT 
~27 



(n|p|n). 



(23) 



(24) 



The same result is obviously obtained starting with the complete set of eigenstates of *. Since we will repeatedly go 
from the eigenstates of * to eigenstates of II or vice versa, let us see how it proceeds 



i 



Because 



one finds 



Tr[p] 

= exp 

(6) 



P*(*|pTI|*) = 



1 



Tr[p] 



p* —iu 2 <* I n x ) (it |p|n 2 ) (n 2 1 *} . 

ZTT Z7T 



i / d A x n* 



= exp 



i / d d x n* 



Tr[p] 



pn x pn 2 

P* — -— -exp 
2tt 2tt p 



i / cfa;(IIi - n 2 )* 



n 2 (ni|p|n 2 ) 



(25) 



(26) 



(27) 



Now one performs the integral over *, which generates ^IT — II 2 ], and after taking the integral either over IT or II 2 
one reproduces the result (|23|) . 
Let us now compute (LT) as 



i r pn 

(H) = - / P* — P^exp 



/ / d 3 x Hip 



n<* + - V j|p|*-- v »). 



(28) 



With the complete sets of momentum eigenstates, the formula (|28| is rewritten as 

■ expjz y'd 3 x[(n 1 -n 2 )* + (i(n 1 + n 2 )-n)^}n(n 1 | / 5|n 2 ). (29) 



if pn pni pn, 
<n ) = - / p* — Vw — - — - 

w Z J 2tt y 2tt 2tt 

The integrals over * and <p produce the deltas S[Hi — n 2 ] and 5[(Ui + n 2 )/2 — II], respectively. After taking the 
trivial integrals over n! and n 2 one shows that 

(n) = (n). (30) 

If one considers an average of a non-local product II(xi)II(x 2 ) . . . II(x„), the result is obviously analogous to Eq. 
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D. C>(*,n) depends on *II 

Let us compute (II(x)*(y)). The computation proceeds in the same way independently whether x = y or x / y. 
So, the arguments are dropped. Using the complete sets of momentum eigenstates, we have 



(II*) = 

Now we evaluate (II*) as 
1 



i r dil vu 2 

— / £>* -exp 

Z 2tt 2tt y 



i / d a x(U! - n 2 )* 



n! + n 2 



*(iii| / s|n 2 ). (3i) 



(n*) = 



Trip] 



D*(*|,3n*|*) 



Tr[p] 



_ T Dili £>n 2 

2?*— — — -exp 
2tt 2tt v 



i / ^(iu -n 2 )* n 2 *(ni|p|n 2 >, (32) 



which is not equal to (II*). 

In case of bosonic fields it was found [6j that (11$) = 5 (A* + $11) that is the operators II and $ have to be 
symmetrized. One can guess that the fermionic operators need to be antisymmertrized. To clarify this point we 
compute (*n) 



(*n) 



1 



Tr[p] 



D*(*|p*n|*) 



1 



Tr[p] 



Z>*(*|f[p*|*). 



(33) 



The minus sign occurs because Tr[AB] = — Tr[.Bj4] when the elements of matrices A and B belong to a Grassmann 
algebra. Using the complete sets of momentum eigenstates, one finds 



(*n) = - 



1 



pit x>n 2 

X>*-^-!--^exp 



Combining the results 



Ti-[p\ J 2tt 2tt 
we prove that 



/ / d 3 x(ni-n 2 )* ni*(ni|/5|n 2 ). 



(n*) = -(*n-n*). 



(34) 



(35) 



The result can be easily generalized to a non-local multiple product of * and II. 

We conclude this section by saying that the Wigner functional ([3]) indeed represents a density in a phase space 
spanned by * and *. 



IV. DIRAC FIELD 



To proceed further a dynamics of fermionic field must be specified. We consider here the four-component Dirac 
spinor *(f,x), which is an operator acting in the Hilbert space of states, with the Lagrangian density 



£ = *(i 7 Al <9 Al -TO)* + £ / , 



(36) 



where * = *^7°, 7^ are the gamma matrices and m is the fermion mass; the interaction term is assumed to either 
describe the coupling to an electromagnetic (^4^) or scalar (*) classical field 



-e*7^*A p , 
-g** * 



with e and g being the coupling constants. The field operators * and * obey the Dirac equations 

25*$. 



As well known, the momentum conjugate to * is 



<9* 



(37) 



(38) 



(39) 
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with ^ = ^ and the field operators obey the anticommutation relations 



{* a (t, x), 9p(t, x')} = 0, {ft a (t, x), ft^t, x')} = o, 
x), xO} = 5^<5 (3) (x - x'), 



where a,(3 = l, 2, 3, 4 are the spinor indices. 
The Hamiltonian is 



d^xU^l 7 • V + im 



and the Hamilton equations, which read 



i 6H m, v- ■ 

* = — ^r- = -7 7 • V + im 

t = - s 4 

59 



ig9 



ie^A^ 
ig<f> 



9, 



9 



(40) 
(41) 

(42) 

(43) 
(44) 



are identical with the Dirac equations 



V. EQUATION OF MOTION 



The equation of motion of Wigner functional is derived from the equation satisfied by the density matrix operator 



i§- t p(t) = [H,p(t)}, 



(45) 



where the Hamiltonian is given by Eq. (|42[) . The field operators are assumed to be in the Schrodinger picture and 
thus the operators are time independent. We split the Hamiltonian as 



where 

H v = - J d 3 xttj a f • V9, 
Eq. (|4"5")) provides 



H — H y + H m + Hi , 



H m ee -im I d 3 x r±7°* 



Hj EE —% \ d 3 xU 



(46) 



9. (47) 



i— W[9,Il;t] = I v + I m + Ij, 



where 



with A = V, m, I. Further on we compute I\r, I m and //. 



h = I Vipexp 



-i \ d 3 xu ip (9 + -<p\[A x ,fi]\Hf- -ip) 



A. Mass term 



(48) 



(49) 



We first analyze the mass term showing some details which will be skipped in the discussion of 7y an d Ii- The 
mass term has the form 



im I d x I Dip exp 



i / d xUip 



(* + i<^l*7o T n/5 + pn 7 *|* - -ip) 



mTr[7°] j d 3 x 5 {3) (x = 0) W[9, H; t], 



(50) 
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where we have used the commutation relation (1411) providing 

n(x) 7 °#(x) = -#(x) 7 Jn(x) + iTr[ 7 °] <5 (3) (x = 0). 



(51) 



The last term in Eq. (|50[) . which is rather pathological because of S^ 3 \k = 0), actually vanishes as Tr[ 7 M ] = 0. It is, 
however, kept here as analogous expressions in Jy and Ij are nonzero. We will show that the pathological term from 
Eq. (1501) will be exactly canceled by another term of this type. It is worth noting that such pathological terms do 
not show up when the equation of motion of Wigner functional of bosonic fields is derived Q . The point is that the 
Hamiltonian of the scalar field is the sum of the terms which depend either on $ or II but not on both. 
Since + j-<p\, \^ — i<p) are eigenstates of ^, we have 



I' m = im / Dip exp 



d 3 3 



i I <rx Hip 
1 



2 rir 1 2 

where the pathological term is excluded from I' m that is 

I' m = I m ~ mTr[ 7 °] f d 3 x8^(x = 0) W[V, U] t]. 
Now, the complete sets of momentum eigenstates are introduced and one gets 

i / d 3 xllip expji / d 3 x (lli - U 



(52) 



(53) 



. , PIT VU 2 

l m — im I Dip — ^— exp 



d 3 3 



2?r 2tt 



* 7o T (ni-n 2 ) + ^i+Ik 



lT , , Hi + n 2 
* + — o — f 



(ni|/5|n 2 ). 



(54) 



Observing that the identity, which is obvious for bosonic fields, holds for the fermionic ones with merely changed 
the sign 



we find 



i— exp 
(5* 



, Dili Dn 2 

I m = -m Vip — — exp 

Z7T Z7T 

5 T S \ 



Ilexp 



d 3 xIW 



(55) 



- i d 3 x U<p 



x / d 3 x(^jQ 
The first term in Eq. (|55|) provides 



exp • 



{z / d 3 ^ [(it -it)* 



Hi + n 2 



-f 



}(ni|/5|n 2 ). 



C = / d 3 x^^w[^,n ; t} 



(56) 



(57) 



but the second one requires further manipulations. 
Using again the identity of the form (|55|) . we get 



I' = -im / Vip 



Dili Dn 2 

2?r 2tt 



M eXP 



i y d 3 a;ll<^ 7 (f-^-exp|z ^ d 3 a 



IT -n 2 * 



ni + n 2 



}(ni|p|n 2 ). (58) 



With the help of partial integration formula, which also holds for the Berezin integrals, we obtain 

C =-mJ d 3 x^nw[^,n ; t}. 



(59) 
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Since 



(7o T ) 



^(x) = _( 70 )^n^(x)— ^-r +Tr[ 70 ]^ 3 )(x = 0), 



where the spinor indices are explicitly written for clarity, we have 



IL„ = m 



JePxItyoj- W[¥, U.]t] -Trfro] J d 3 xS^(x = 0) W[%U;t]. 



(60) 



(61) 



One sees that the pathological term encountered in Eq. (|50p is exactly compensated by the term found in Eq. (foT 
Therefore, combining the results (|57[ 1611) . the final expression of I m reads 



m / d x 



t S „ 5 



W[*,Il;t]. 



(62) 



B. Interaction term 



The computation of the interaction term J/ proceeds in exactly the same way as that of the mass term I m but 771.7° 
must be replaced by either e 7 °7 / M Al or (77°$. Therefore, we have 



h 



d 3 s 



e 7 °7 A1 ^ 
37°$ J M> 



ff 7°$ J 5U 



W[%U]t]. 



(63) 



C. Gradient term 

It is more tedious to compute the gradient term Jy than the mass one I m but conceptually there are no important 
differences. Finally, one finds 



i / d x 



W[*,Il;t], 



(64) 



where the relative minus sign of the two terms results from the additional, when compared to analogous terms in J m , 
partial integration with respect to x. 



D. Equation of motion 



Summing up the terms Jy, I m an d Ii, the equation of motion of Wigner functional (|48[) gets the form 



at 



W[9,U;t] 



* 7 7 • V + 77777° + i 



d 6 x 



IT 7°7- V -77777° 



e 7 7 - 



57° $ J J 611 

Eqs. (|43l |44|) show that the equation of motion can be rewritten as 

SH 5 5H 5 



_6_ 



W[^,H;t] = 0. 



d_ 

di 

where the 'classical' Hamiltonian equals 



d 6 x 



SU a S^a S^ a SU a 



H 



d 6 x 



n 7 ° 7 • v + 7777 + 7 



g<f> 



= / d 3 x 



W[*,IL;t]=0, 



n- ll (7-v-/ ( ii-i[ e ^" 



(65) 



(66) 



(67) 



We note that the differentiation of the 'classical' Hamiltonian with respect the 'classical' fields proceeds somewhat 
differently than the differentiation of the Hamiltonian operator with respect to the field operators, as the 'classical' 
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fields belong to the Grassmann algebra. Therefore, the signs of 5H/5U and SH/5H coincide but those of 6H /Sty and 
differ. 

The equation of motion (|66j) has a structure of the classical Liouville's equation. This is an important difference 
when compared to the Wigner functional of bosonic fields where the Liouville's equation is obtained only when the 
fields are either free or the quantum corrections to the interaction are neglected [6J. The fermionic field instead 
remains always classical except the effects related to the Pauli principle. This is not surprising as the Dirac equation 
is linear in the fermion field. 

Once the equation of motion (|66l) is of the Liouville's form, one immediately writes down its 'elementary' solution 

as 



Wei[*,n ; t] = $[*(x) - * c (t,x)] <5[n(x) -n c (t,x)], (68) 

where v I / c (t, x) and U c (t,x) is the solution of 'classical' equations of motion 

t , x SH • , , 6H 

*^ = m^y n ^ x) = M^)> (69) 

which coincide with the Dirac equation and its conjugate. We note that the sign in the second equation is different 
than usual, as the differentiation is performed in the Grassmann algebra. The solution (|68| is called 'elementary' 
because it corresponds to a specific solution *f? c (t,x), n c (i,x) of the equations of motion. A general solution of 
Eq. (j6"6]l is a superposition of the elementary solutions (1551) with various \Ij c (i, x), H c (t, x). One typically considers a 
solution of Liouville's equation which includes a distribution of initial values of the elementary solution. 



VI. FINAL REMARKS 



Quantum mechanical problems are usually solved in terms of wave functions determined by the Schrddinger equa- 
tion. The approach using the Wigner function, which obeys the transport-like equation of motion, offers an alternative 
technique. However, it has taken quite some time to develop appropriate methods 043- The counterpart of the for- 
malism employing the Wigner functional is the Schrddinger functional equation satisfied by the wave function which 
functionally depends on a field. The approach, which was first introduced for the bosonic fields [|[ and then extended 
to fermionic ones [ToL Hlj |. has various applications. Nonperturbative studies of gauge fields and quantum gravity 
can be enlisted here, see [l2|, EH for exemplary recent developments. Hopefully, methods referring to the Wigner 
functionals will be also useful. 
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